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ABSTRACT 

We have obtained high spatial resolution Keck OSIRIS integral field spectroscopy of four z ~ 1.5 
ultra- luminous infrared galaxies that exhibit broad Ha emission lines indicative of strong AGN activity. 
The observations were made with the Keck laser guide star adaptive optics system giving a spatial 
resolution of 0.1" or < 1 kpc at these redshifts. These high spatial resolution observations help to 
spatially separate the extended narrowdine regions — possibly powered by star formation — from the 
nuclear regions, which may be powered by both star formation and AGN activity. There is no evidence 
for extended, rotating gas disks in these four galaxies. Assuming dust correction factors as high as 
A(Ha) = 4.8 mag, the observations suggest lower limits on the black hole masses of (1 — 9) x 10 8 
Mq, and star formation rates < 100 M Q yr _1 . The black hole masses and star formation rates of the 
sample galaxies appear low in comparison to other high-z galaxies with similar host luminosities. We 
explore possible explanations for these observations including, host galaxy fading, black hole growth, 
and the shut down of star formation. 

Subject headings: instrumentation: adaptive optics — galaxies: active — galaxies: high-redshift - 
black hole physics 



1. INTRODUCTION 

Spitzer 24 /im imaging in the 9 square degree NOAO 
Deep Wide Field Survey CNDWFS: IJannuzi fc DevlfTool 
of Bootes has revealed large samples of high red- 
shift ultraduminous infrared galaxies (ULIRGs; Liji — 
10 12 ~ 13 ). A simple optical-to-IR color cut of R - [24] > 
14 Vega mag, i.e. 7^(24 ^m)/ f v {R) > 1000, selects for 
an extre me class of ULIR Gs at redshifts of < z > — 1.99, 
a z = 0.5 (|Dev et al.ll2008D . ULIRGs selected this way are 
known as Dust Obscured Galaxies (DO Gs), and are typ i- 
cally redder than low redshift ULIRGs (TDev et al.l f2008). 
The extreme luminosities and colors of DOGs suggest the 
presence of AGN accretion and/or intense star- formation 
heavily obscured by dust at rest-frame optical and UV 
wavelengths. 

The rest-frame optical-to-mid IR (MIR) spectral en- 
ergy distributions (SEDs) of the DOGs show two classes: 
(1) power law sources with SEDs that rise smoothly into 
the MIR; and (2) so called "Bump" sources, with a rest- 
frame 1.6 /im peak in their SED thought to be produced 
by the photospheres of cool stars. The power-law DOGs 
have been shown to exhibit strong AGN characteristics 
including broad-Ha (jBrand et al.l 120071 ) , and a lack of 
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PAH emission in MIR spectra from Spitzer ([Houck et al.l 
I2005D . Meanwhile MIR spectra of bu mp DOGs tend to 
show PAH emission (jDesai et alll2008[ ) , a good indication 
of vigorous ongoing star formation. These two classes of 
DOGs appear to differentiate with luminosity, with the 
more luminous sources (e.g. F 2 4 > 0.75 mJy) having 
a higher frac tion of power-law SEDs likely to harbor ob - 
scured AGN (|Magliocchetti et al.|[2007l: iDev et aLll2008fl . 

The space densities and clustering strength 
of the DOGs are simil ar to sub-mm g alaxies 
CSMGs: iBlairi et al l [2001 iMagliocchetti et all [20071 
iBrodwin et all 120081 ). which are thought to have merger 
driven star formation rates as high as 1000 Mq yr _1 . 
The clustering strengths are also similar to massive 
elliptical galaxies in the local universe, suggesting a 
possible connection between these three different galaxy 
classes. 

Detailed merger simulations of massive gas rich galax- 
ies have been shown to repro duce the selection crite - 
ria of both SMGs and DOGs ([Narayanan et all 12009a! ) . 
In these simulations, the SMGs are typically classi- 
fied as bump sources, and the SMG phase is fol- 
lowed by a hotter, post-merger, AGN dominated phase 
where the galaxy is tran sformed into a power-law DOG 
([Narayanan et all l2009a|) . Currently, mergers are the 
only model shown to produce the extreme colors and 
luminosities of DOGs. 

Because of their potential link to gas rich merg- 
ers, it was surprising to find that the rest-frame opti- 
cal morphologies of the more luminous DOGs are typ- 
ically smooth with little sign of an ongoing merger 
(as evidenced by doub le nuclei. IMelbourne et all [2008b: 
iBussmann et all 120091 : IMelbourne et all 120091 ) In fact, 
the bulk of DOGs with high spatial resolution imaging 
(from HST or Keck adaptive optics) show disk-like or 
elliptical like morphologies, rather than double nuclei or 
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tidal tails. However, their physical half-light sizes are 
small in comparison with LIRGs at z ~ 1, which also 
favor disk-like profiles; the more luminous DOGs have 
typical half-light sizes that are a fact or of 2 (or more) 
smalle r than sample of z ~ 1 LIRGs (Melb ourne et al] 
2008a, h). The unusual sizes of the DOGs are consistent 
with the post-merg er products of gas rich me rgers in the 
local universe (e.g. Rothbcrg & Joseph 2006). 

If DOGs have undergone recent merger activity, evi- 
dence for the merger should be imprinted on the kine- 
matics of the gas in the galaxy. While they exhibit expo- 
nential disk profiles, the gas in the DOGs may not nec- 
essarily be undergoing ordered rotation. With its high 
spatial and spectral resolution, the Keck OSIRIS integral 
field spectrograph (IFS) is a ideal tool for measuring gas 
kinematics in high-z galaxies. OSIRIS is an adaptive op- 
tics fed instrument with a diffraction limited angular res- 
olution of 0.05" at A-band (2.2 /mi). OSIRIS produces 
a spectrum at every spatial resolution element over a 
field of view of 2 — 3 arcsecs. With a spectral resolution 
of R ~ 3000, OSIRIS can resolve kinematic signatures 
with FWHM of 100 km s -1 . Targeting the Ha emis- 
sion line redshifted into the near-infrared (NIR), OSIRIS 
has been used to de termine the kinematics of z = 1 — 3 
galaxy samples (e g. iWright et aLll2007t lLaw et al.ll200l 
IWright et al"]l2009h . 

Because of its high spatial resolution, OSIRIS also 
offers the possibility of spatially differentiating the ex- 
tended star forming region from the central regions which 
may contain both star formation and AGN activity. This 
capability allows for much cleaner determinations of the 
relative strengths of these two power sources, compared 
with seeing limited observations. 

This paper presents OSIRIS observations of four DOGs 
at z > 1. Each DOG was targeted at the wavelength of 
the Ha line redshifted into the NIR. We quantify the flux 
contributions from broad and narrow Ha emission-line 
regions, and use the data to place constraints on super- 
massive black hole masses, star formation rates (SFRs), 
and metallicities of the DOGs. We also examine the kine- 
matics of each system. 

This paper is organized as follows. The sample selec- 
tion and OSIRIS observations are described in Section 2. 
Section 3 details the Ha line measurements. Properties 
of the sample derived from spectra, including the black 
hole masses and star formation rates, are presented in 
Section 4. The results are discussed in Section 5, and we 
compare the DOGs to other high-z and local AGN hosts. 
Conclusions are provided in Section 6. Throughout we 
assume a canonical A Cold Dark Matter Universe with 
Cl M = 0.3, fl A = 0.7, and H Q = 70 km s" 1 Mpc^ 1 . 

2. OSIRIS INTEGRAL FIELD SPECTROSCOPY 

We used the Keck laser guide star adaptive optics 
(LGSAO) system and the OSIRIS IFS to observe four 
z > 1 ULIRGs. These galaxies were selected from op- 
tical and MIR photometry of the NDWFS of Bootes to 
have R - [24] > 14 Vega mag, and are thus DOGs. The 
sample selection also required that each DOG have a 
strong Ha detection in seeing limited NIR spectroscopy. 
These spectra provided not only a redshift, but also an 
Ha flux. We selected those DOGs with the strongest Ha 
detections that were sufficiently near to an AO tip-tilt 
guide star, to provide excellent AO performance. 



The requirement of large Ha fluxes (Fhu > 1 X 10 16 
ergs cm~ 2 s _1 ) was made to ensure that the OSIRIS ob- 
servations were successful. This flux limit translates to 
an Ha luminosity of 5 x 10 41 ergs s" 1 at z = 1 and 
3 x 10 42 ergs s -1 at Z = 2. This flux requirement also 
introduces a significant selection bias. First, DOGs at 
rcdshifts where Ha falls between the H and A"-bands 
(i.e. z ~ 1.9) will not be in the sample. Many of the 
Bump DOGs lie in this redshift range because it places 
the 8 /im PAH feature into the Spitzer 24 /im band, red- 
dening the the optical-to-MIR color. Thus all of the 
DOGs in our sample have power-law SEDs (Figure [lj 
and are likely to be powered by AGN. In addition, by 
selecting sources with the highest Ha fluxes, the sample 
skews to somewhat lower redshift than the typical DOG 
(which has < z >^= 1.9 ±0.5). The final sample also 
skews to higher 24 /im flux with 3 of the 4 DOGs having 
F„(24) > 1 mJy 

2.1. The Observations 

OSIRIS is an AO-fed instrument. The Keck LGSAO 
system uses a deformable mirror to correct for atmo- 
spheric distortions to the wavefront, and thus recov- 
ers the diffraction limited resolution of the telescope 
(~ 0.05" at 2.2 /im). The system tracks atmospheric 
distortions with observations of a reference star within 
50" of the science target (for tip-tilt correction), and a 
sodium laser guide star propagated from the telescope to 
the location of the target (for higher order corrections). 
Table [2] summarizes the conditions during OSIRIS ob- 
servations, including magnitudes and separations of the 
tip-tilt guide stars. 

Observations taken in 2008 were made under photo- 
metric conditions with the instrument at the optimum 
well controlled temperature. Unfortunately, during 2009, 
OSIRIS began to heat up, 5 degrees (K) above the op- 
timal operating temperature. The higher temperature 
produced noisier data and resulted in several issues that 
made 2009 data more difficult to reduce: the dark current 
was about 40% higher compared to 2008, and warmer op- 
tics made extracting high signal-to-noise spectra more 
difficult. The observing conditions were clear during 
2009 observing runs. 

In all cases, we observed with the 0.05" plate scale us- 
ing narrow spectral filters centered on Ha, thereby allow- 
ing the largest possible field of view for the chosen plate 
scale. The typical field of view was 2" x 3" or roughly 
20 x 30 kpc at these redshifts. 

Because of the small field of view of OSIRIS, we first 
centered the AO tip-tilt guide star onto the OSIRIS 
frame. We then manually offset to the target using off- 
sets derived from NDWFS imaging. This method pro- 
vided good acquisition of the galaxy on the OSIRIS field 
of view. 

Individual exposure times were 15 minutes with total 
on-source exposures ranging from 1 to 3 hours. Small 
dithers were applied after each exposure. However, the 
galaxy was always within the field of view of the instru- 
ment, and the galaxies were sufficiently small that no 
additional sky frames were needed. 

2.2. Data Reduction 

The data were reduced with the OSIRIS reduction 
pipeline. For each set of observations, a sky frame was 
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TABLE 1 

DOG Properties Summary 



Galaxy 


Object Name 


redshift 


M24) 


R 








[mjy] 


[Vega mag] 


NDWFS.J143027.2+344008 


DOG1 


1.370 


1.169 


24.76 


NDWFS.J143335.6+354243 


DOG2 


1.300 


5.577 


21.85 


NDWFS.J143400.3+335714 


DOG3 


1.684 


1.754 


23.89 


NDWFS.J143424.5+334542 


DOG4 


2.260 


0.860 


25.60 



TABLE 2 

Keck Observation Run Summary 



Object Name UT Date 



Exptime Filter 

[s] Center [fJ.m] Width [/jm] 



R 



Guide Star 
[mag] offset 



Observing 
Conditions 



DOG1 
DOG2 



DOG3 
DOG4 



Jun 04, 2008 
May 07, 2009 
May 26 
May 27 
May 07 
May 27 



2009 
2009 
2009 
2009 



4500 
2700 
3600 
3600 
7200 
10800 



1.571 
1.504 



1.765 
2.175 



0.078 
0.075 



0.087 
0.108 



16.3 
15.1 



16.0 
13.4 



31 
47 



39 
32 



some cirrus 
clear 
clear 
clear 
clear 
clear 



created by median combining dithered science frames. 
Each frame was sky subtracted, and cleaned of cosmic 
rays. A 3D data cube was then extracted using rectifi- 
cation matrices supplied by Keck Observatory. Frames 
were mosaicked together to produce a final 3-D data 
cube. 

As a result of running at higher temperature in 2009, 
optical elements within OSIRIS expanded, slightly alter- 
ing the optical path. Therefore the canonical rectification 
matrices used to convert the 2-D raw frames of overlap- 
ping spectra into 3-D data cubes were no longer valid. 
While new rectification matrices were produced they did 
not work as cleanly as the original, and when applied 
to the data they created some pixels of unrealistically 
high and low count levels (e.g. ±10 u counts). In order 
to remove these glitches we used a median combine (as 
opposed to a mean) when making the final data cubes. 

After each science target, an A0V star was observed 
with an identical instrument setup. These observations 
were used to correct the spectra for Telluric absorption 
and to flux calibrate the images. After the basic reduc- 
tion of the standard star observations, a ID spectrum was 
extracted of the star. Hydrogen absorption lines were 
modeled out using the OSIRIS pipeline Telluric Extract 
routine. The stellar spectra were divided by a black-body 
spectrum with a temperature of an A0V star. The re- 
sulting spectrum was then divided into the science data 
to correct for Telluric absorption. 

The spectra of the standard stars were also used to 
apply a rough flux calibration to the data. A total flux 
for each star was measured from the OSIRIS data cube 
and compared to the expected flux of the star across 
the OSIRIS filte r given the J, H, and K band fluxes 
(jElias et al.lll982h . To determine the expected flux of the 
standard star in the OSIRIS filter we assumed a black- 
body spectrum across the filter band-pass, and assumed a 
flat filter function. To check this calibration, we compare 
the flux in Ha for DOG1 (See Table 1) to that found in 
iBrand et~a l. (2007) from long slit spectroscopy. Our Ha 
measurement is 60% larger than the Brand et al. value, 



which is not surprising as they suggest that there data 
could suffer slit losses on that order from pointing issues. 
Because OSIRIS is effectively imaging in Ha, slit losses 
are minimal. 

3. SPECTRAL LINE MEASUREMENTS 

Figure[2]shows the Ha maps of the four DOGs observed 
by our program. All show point-like Ha morphologies. 
We center a circular extraction region on the peak of 
the Ha flux, and extract a region with radius of 0.15", 
enclosing the diffraction limited core of the PSF which 
is measured to have a FWHM of 0.1". The spectra of 
these point-like objects (Figure [2]) show broad Ha emis- 
sion. We first fit these broad lines with a single Gaussian 
profile (red curve in Figure[2]). The resulting fits have line 
widths in excess of 2000 km s _1 indicative of AGN ac- 
tivity. While the spectra are noisy, we also attempt a fit 
that includes narrow-line Ha and [Nil A6548,6583]. The 
total fluxes and line widths of the broad Ha line do not 
change significantly with the inclusion of these narrow- 
lines in the fit. Fluxes and line- widths of the broad Ha 
emission are given in Table [3] 

We also attempt to identify the spatial extent of the 
narrow-line emission in each system. To accomplish this, 
we first subtract off the broad- line AGN component. We 
create a point-spread-function (PSF) image by collaps- 
ing the data cube in the spectral region that encompasses 
the wings of the broad-line Ha spectrum. At each wave- 
length, we scale the peak of the PSF image to the flux 
of the modeled broad Ha line and subtract it from the 
data cube. The result of this process is a 3-D data cube 
of the narrow Ha flux with the broad-line Ha removed. 

The narrow Ha images of each galaxy are shown in Fig- 
ure [3] These images have been smoothed by a Gaussian 
kernel with FWHM of 2 pixels to enhance the detection 
of any prominent sources. In general, these images do 
not show regular morphologies such as those expected 
if star formation were smoothly distributed within well- 
behaved disk structures. Rather, these regions appear as 
individual clumps. Except for DOG2, these narrow-line 
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D0G1 NDWFS_J1 43027.2 + 344008 



DOG3 NDWFS_J1 43400.3 + 33571 4 
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Fig. 1. — Optical through MIR SEDs of the four DOGs observed with OSIRIS. These SEDs rise steadily into the MIR, and do not contain 
features such as a 1.6 /im stellar bump, or a Balmer break. High redshift ULIRGs with SEDs like these have been generally shown to 
contain dust obscured AGN. 



regions are separated from the nuclear region which hosts 
the AGN (indicated by green contour lines). 

Shown on the images are the regions used to measure 
the narrow Ha line widths and strengths. Some regions 
that appear bright in the narrow Ha images were not in- 
cluded in the creation of the final spectrum. Individually, 
the spectra of these regions showed no indication of con- 
taining emission line flux. Rather, these regions appear 
to be noise fluctuations in the collapsed 2D images. 

The spectra of the narrow-line regions are shown in 
Figure [3] Evidence for narrow Ha emission is found in 
each DOG. We fit the spectra for Ha and [Nil] lines, 
using Gaussian profiles (red line in Figure). FWHMs of 
the narrow emission lines vary from ~ 200 — 500 km s -1 , 
with typical fluxes from (0.5 — 2.5) x 10~ 17 ergs cm -2 
s . The velocity offsets between the narrow and broad 
Ha lines in each galaxy range from —150 to 500 km s _1 . 

Three of the DOGs show evidence for [Nil] emission. 
DOG1 and DOG3 have [Nil] /[Ha] line ratios of - 0.3. In 
contrast, DOG2 has [Nil]/ [Ha] = 0.74. Such a high ratio 
likely indicates that the production of [Nil] in this galaxy 
is be ing enhanced by AGN activity (jSwinbank et al.1 
12004J ) , a reasonable assumption given that the emission is 
located spatially atop the nuclear broad Ha flux. Shocks 
in the nuclear gas could also result in high [Nil] /[Ha] 
ratios. The narrow-line measurements are summarized 
in Table H 

4. PROPERTIES OF THE SAMPLE 

4.1. Dust Extinction 

In nebular theory, the Ha/H/3 ratio is fixed by 
the radiation field. For case B recombination, in a 
density-bounded HII region powered by star formation, 



Ha/Hff= 2.86. Fo r AGN-ionized regions, Ha/H/3= 3.1 
(jOsterbrockl Il989f ). Ha/H/3 ratios larger than these 
canonical values indicate the presence of dust, which will 
preferentially attenuate the bluer H/3 line compared with 
the Ha line. Thus, a measurement of the Ha/H/3 line ra- 
tio can be used to estimate the reddening in a system. 

The OSIRIS observations of the DOGs in our sample 
were only made in a small spectral window (~ 100 nm 
wide) about Ha and therefore do not contain the infor- 
mation necessary to produce an estimate of dust extinc- 
tion. However, 3 DOGs in the sample were observed pre- 
viously with NIR long-slit spectrographs covering both 
Ha and H/3 lines. While H/3 was not strongly detected 
in any of these observations, limits on the line strengths 
were obtained. 

DOG1 was observed with the Gemini N IRI spectro- 
xaph and the results were published in iBrand et al.1 
20071 object J143027.1+344007). In this system, the 
broad-line Balmer decrement Ha/ H/3 > 22.5, which, as - 
suming a Milky Way dust curve (ICardelli et al.lll98l . 
translates to a reddening, E(B — V) > 2.09 mag, or 
A(Ha) > 4.86 mag. Thus the extinction corrected Ha 
emission in this system is likely to be a factor of 70 larger 
(or more) than measured. 

DOG3 and DOG2 were observed by Palomar Triple- 
Spec (Melbourne et al. in preparation). Unfortunately, 
the limits on the Ha/H/3 ratios in the Palomar data are 
less stringent than the NIRI data, with Ha/H/3 > 6.5 and 
7.6 respectively. These are equivalent to A(Ha) > 1.84 
mag and 2.22 mag, or roughly a factor of 5 - 10 in flux. 
However, A(Ha) could be significantly larger then this 
lower limit. 

For DOG4 we do not have a measurement of H/3, and 
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Fig. 2. — Ha maps and ID spectra of the sample ULIRGs. Images are 1 arcsec on side (~ 10 kpc at these redshifts). All four DOGs 
show a point-like morphology in the Ho flux distribution. This point-like flux is dominated by a broad Ho line, with FWHMs > 2000 km 
s —1 , suggesting AGN activity within the DOGs. 

therefore do not have a measure of the dust extinction. 
We choose to adopt the minimum dust extinction mea- 
sured from the other 3 DOGs (A(Ha) > 1.84 mag) as a 
lower limit on the extinction for DOG4. 



4.2. Black Hole Masses 

Under certain assumptions the OSIRIS data cubes can 
be used to estimate the black hole masses in our sam- 
ple galaxies. To make the estimate we require two mea- 
surements: (1) the Ha line width, and (2) a continuum 



flux from the AGN point source. The hrst quantity we 
measure directly from the spectra. Unfortunately, the 
second quantity is difficult to measure directly from the 
OSIRIS data. While continuum is seen in the spectrum 
of the central region of the galaxy, it is not clear how 
much of this continuum is produced by the AGN and how 
much is produced by the galaxy. However. IGreene fc Hoi 
(2005) show that there is an empirical relationship be- 
tween broad-line Ha flux and the continuum flux of the 
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DOG1 




Fig. 3. — Same as Figure [2] only now the broad line component has been subtracted from the image, revealing the narrow-line Ha flux. 
These images have been smoothed by a Gaussian kernel of 2 pixels width to bring out the contrast. Each DOG appears to contain one or 
more narrow-line regions often spatially offset from the broad-line region (contours). Cyan boxes show the regions included in the extraction 
of the narrow-line flux. These regions tend to be clumpy rather than extended (as might be expected in a relaxed disk). The top axis 
shows the velocity offset from the broad-line AGN component. The spectra also show evidence for narrow [Nil] lines. 
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AGN, 



L5100 = 1.23 x 10 7 ■ (L Ha 



0.864 



(1) 



We use the reddening-corrected broad-line Ha luminosi- 
ties to determine the AGN continuum luminosity at 5100 
A, L5100. With the Ha line width and continuum flux, 
we calcula te the black hole m ass from the empirical cal- 
ibration of iPeng et all (|2006al ). 



L5100 



1 x 10 44 [ergs/s] 



0.5!) 



line-width 
1000 [km/s] 



2.0(i 



M Q (2) 



The estimated black hole masses are summarized in Ta- 
ble [3] they range from (1 - 9) x 10 8 M Q . DOG3 hosts 
the most massive BH, not surprising given the width of 
the broad Ha line in that system. The intrinsic scatter 
in this empirically derived relation (E quation [21 is about 
0.3 — 0.4 dex (e.g. Vestergaard 2002). However, system- 
atic uncertainties in the dust correction could skew these 
results to lower BH masses and all reported BH masses 
are lower limits. 

4.3. Rest-Frame Optical Galaxy Stellar Luminosity 

To better compare the DOGs with other galaxy sam- 
ples (Section 5) we attempt to measure the stellar con- 
tribution to the galaxy luminosity, independent from the 
AGN contribution. Unfortunately, because the spatially 
extended galaxy continuum is not well detected in the 
OSIRIS frames, we can not use the OSIRIS data to con- 
strain the relative contribution of AGN and star-light 
to the total luminosities. However, we get some hints 
from high spatial resolution HST and AO imaging of 
larger DOG samples, which show that the AGN (cen- 
tral point source) typically contributes only ~10% of 
the galaxy light at NIR (rest-frame optical) w avelengths 
(Me lbourne et al.ll2009t IBussmann et al.ll2009fl . 

While this ratio could prove to be a useful rule of 
thumb for DOGs, we know that it does not hold univer- 
sally. In fact, we have a high spatial resolution AO im- 
age of DOG 1, which is actually point-source dominated 
^Melbourne et al.ll2009D . suggesting that the AGN could 
be contributing a much higher fraction of the light in 
this system. This example suggests that in the absence 
of imaging we should attempt some other method for 
determining the AGN contribution to the galaxy lumi- 
nosity. 

Fortunately, the OSIRIS data can be used to mea- 
sure the AGN continuum contribution based solely on 
the Ha flux. As quoted from Greene & Ho, the flux 
from Ha correlates with the AGN continuum emission at 
5100A(L5100). With a measure of the AGN continuum 
in hand, we convert to total AGN flux in a given pass- 
band assuming a standard QSO spectrum. Subtracting 
the AGN flux from the total observed flux (at rest-frame 
optical wavelengths), we derive the stellar luminosity of 
each galaxy which we convert into a rest-frame R-band 
absolute magnitude assuming a power-law spectral en- 
ergy distribution for the galaxy light. As expected from 
its high spatial resolution AO image, the flux of DOG1 
is dominated by AGN light; the AGN contributes 95% of 



the total. However, the other DOGs show significantly 
less contribution from AGN light to the total luminosity, 
~ 10%, which is typical for DOGs based on results of 
image decomposition. 

Un-corrected for reddening, the DOGs in our sample 
have Mr = —20 to —24. However, DOGs are among the 
most dust obscured galaxies known, with colors redder 
than the typical low-z ULIRGs. This suggests that even 
the galaxy light may be heavily attenuated by dust. If 
we apply the dust correction from the broad-line regions 
to the entire galaxy (likely an upper limit on the actual 
dust obscuration for the galaxy) then the DOGs are very 
luminous, M R = -24 to -26.5. 

4.4. Star Formation Rates 

After subtracting the broad- line Ha component, each 
DOG shows additional narrow-line Ha flux. If we as- 
sign all of the observed narrow Ha flux to star formation 
we can estimate the star formatio n rate (SFR) o f each 
galaxy. We u se the prescriptio n in iWright et al.l (|2010l ) 
adopted from iKennicuttl (|1998f) . 



SFR [M /yr] = L(Ha)/2.5 x 10 
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(3) 



Uncorrected for reddening, we derive typical star SFRs 
for the DOGs of only 0.5 - 2.0 M Q yr" 1 (Table H}. 

As with the broad-line fluxes, these narrow fluxes are 
likely to be heavily dust extincted. Because they are 
spatially removed from the center of each galaxy, the dust 
obscuration may be different than for the AGN. However, 
we can use the AGN obscuration as a rough estimate. 
Doing so gives SFRs 1-2 orders of magnitude larger than 
measured. Even with these correction factors the star 
formation rates that are measured are small for ULIRGs, 
< 100 Mq yr _1 . In addition, some of the narrow-line flux 
could be produced by the AGN, meaning that the SFR 
could be significantly lower than these limits. 

In addition, we have shown that we are sensitive to 
very small star formation rates. Clearly the surface den- 
sity of star formation in any extended component must 
be small, e.g. < 0.5 M0yr _1 kpc -2 (before correcting for 
dust attenuation), or our observations would have de- 
tected it. 

4.5. Metallicity 

Three galaxies in our sample show evidence for an [Nil 
A65831 collisionally e xcited emission line. As shown in 
Melbourne k Salzerl (|2002f) . the [Nil] /Ha line ratio can 
be used as a rough proxy for metallicity in star forming 
galaxies. Assuming the [Nil] emission in these systems 
is entirely the product of star formation we can estimate 
the metallicity of these systems using, 

12 + log(0/H) = 9.26 + 1.23 ■ log([NII]/Ha) + 

0.204- [log([NII]/Ha)] 2 (4) 

For the three galaxies for which [Nil] is observed 
(DOGS 1, 2, and 3), the estimated metallicity is roughly 
solar or larger. However, for DOG2, the [Nil] emission 
lies atop the nucleus of the galaxy (Figure 2) and is likely 
boosted by AGN heating. In this system, [NII]/Ha> 0.7, 
another indicat ion that the [Nil] li ne may be enhanced by 
AGN activity (|Brand et al.ll2007fh Therefore, Equation 
4 is expected to over-predict the metallicity in DOG2. 
The metallicities are given in Table 4. 



TABLE 3 
Broad-Line Ha Properties 



Object Name 


Obs. Ha flux 
[10 — 16 ergs cm -2 s _1 ] 


Ha width 
[km s ] 


Ha Lum. a 
[10 42 ergs s" 1 ] 


A[Ha] 
[mags] 


Black Hole Mass b 
[10 6 M Q ] 


Host Galaxy 
M fl c 


DOG1 


5.98 ±0.25 


2068 ± 22 


6.77 


>4.86 d 


> 490 


-20.06 


DOG2 


1.97 ±0.09 


2288 ± 40 


1.98 


>2.22 e 


> 103 


-23.85 


DOG3 


4.97 ±0.15 


6757 ± 96 


2.60 


>1.84 e 


> 878 


-24.73 


DOG4 


1.36 ±0.20 


2358 ± 221 


5.34 


>1.84 f 


> 138 


-24.30 



a Uncorrected for dust attenuation 
k Corrected for dust attenuation 

c Corrected for AGN contamination, but no dust correction 

d |Brand etTal! d2007l) 

c Melbourne et al. in Preparation 

^ No measurement of H/3 available so used the minimum correction from the others in the sample. 



TABLE 4 
Narrow-Line Ha Properties 



Object Name Ha flux Ha width Ha Lum b SFR b [Nil] /Ha Metallicity c 

[10" 17 ergs cm" 2 s" 1 ] [km s" 1 ] [10 41 ergs s" 1 ] [M Q yr" 1 ] [12+log(0/H)] 



DOG1 2.5 ± 0.2 533 ± 47 2.80 1.3 ± 0.1 0.27 ± 0.03 8.6 

DOG2 a 0.50 ±0.12 470 ± 48 0.50 0.22 ± 0.02 a 0.74 ± 0.24 9.1 a 

DOG3 1.3 ±0.1 210 ±56 2.41 1.1 ±0.1 0.29 ±0.04 8.7 

DOG4 1.4 ± 0.1 152 ± 80 5.53 2.5 ± 0.1 



a The narrow lines observed for TDOG_5 arise from the same region as the broad lines, and the line ratios indicate that these 
lines arc produced by the AGN rather than star formation. As a result the calculated star formation rates and mctallicitics 
for TDOG_5 arc likely incorrect. 

Uncorrected for dust attenuation. 
c Derived from the [NII]/Haratio IMelbourne &: Salzcr 2002), with an uncertainty of roughly 0.3 dcx. 



4.6. Kinematics 

Ideally, we would have detected Ha across an extended 
star forming disk in each galaxy. We could then have 
then used the observed kinematics to model the mass 
distributions within the DOGs. However, due to poor 
S/N and/or the nature of the DOG kinematics, none 
showed evidence for well ordered rotation in an extended 
gas disk. 

Despite these issues, there do appear to be narrow Ha 
emitting regions detected in each galaxy. In DOGs 1 and 
2, the narrow Ha flux arrises from a single source near 
the nucleus, each only 1 — 2 kpc in diameter. The line- 
widths of the narrow Ha emission in these two DOGs 
are resolved at 530 and 470 km s _1 respectively. If the 
ionized gas is tracing a virialized mass distribution, the 
masses within these knots are extremely large, M > 10 11 
Mq. More likely, the gas is not virialized. The line- 
widths could be enhanced by the nearby AGN, or from 
kinematically disturbed gas, or from large bulk motions 
such as from expanding bubbles from a wind. In the 
case of DOG2, the narrow emission resides exactly atop 
the location of the broad-line region, suggesting some 
fraction of the narrow-line flux may be from the AGN. 
The high [Nil] /Ha line ratio of this region also suggests 
AGN contamination. For DOG1, the [Nil] /Ha line ratios 
are more consistent with star formation. 

The other two galaxies, DOG3 and DOG4, appear to 
show extended narrow Ha emission in several distinct 
knots with separations as far as 5 kpc. Each of the knots 
has a very small Ha line width of ~ 150 — 200 km/s, and 
there are no large velocity gradients across these multiple 
emitting regions. 



In each case, the central wavelength of the narrow- 
line component is offset from the broad-line component. 
These kinematics offsets range from 150 km s _1 , in DOGs 
2 and 3, to 350 and 500 km s" 1 , in DOGs 1 and 4 re- 
spectively. 

5. DISCUSSION 

In this set of four DOGs, the broad Ha lines associ- 
ated with AGN have fluxes 10-20 times larger than the 
detected narrow lines produced in star forming regions. 
The relative strengths of these lines suggest that the bulk 
of the dust heating results from AGN activity rather than 
star formation. Based on the widths of the broad lines, 
and the continuum flux of the nuclear region, we esti- 
mate lower limits on the black hole masses for the four 
DOGs in our sample, in the range of (1 — 9) x 10 s M . 

5.1. Comparing the BH Masses of DOGs to BH Masses 
of Other z ~ 2 Galaxies 

Figure [4] compares the BH masses of the the DOGs 
(blue limits) wi th the BH masses in other z ~ 2 and 
z re galaxies (IPeng et al.M 2006bt IRidgwav et "all 120011: 
IKukula et ail 120011) . The z ~ 2 galaxies are a mix of 
radio-loud and radio-quiet quasars. Their black holes 
masses were measured with a technique similar to the 
one used in this paper, except that the emission lines 
used w ere from the rest-fr ame UV, C IV A1549 and Mg II 
A2798 (|Peng et al.ll2006al ). The z ~ 2 host galaxy mag- 
nitudes of the quasars were measured from high spa- 
tial resolution NIR imaging from HST. The compari- 
son sample of interest at z = are elliptical galax- 
ies with BH masses measured from high spatial reso- 
lution HST spectroscopy of the circum-nuclear gas and 
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Fig. 4, — Left: Estimated blac k hole mass vs. host luminosity for the DOGS (blue limits), z ^ 2 AGN (green shading Peng ct al. 2006a; 
Ridgway ct al. 2001; Kukula ct al. 2001) and local galaxies (grey shading). Also shown are X-ray derived BH masses for two of the DOGs 
(red points). In both cases the X-ray results return BH masses that are larger than the OSIRIS derived lower limits. The best fit local 
relation is shown as a solid black line, while the best fit z ~ 2 relation is shown as a dotted line; the dashed lines show different SED 
assumptions used to derive the host galaxy luminosities. Both the BH mass estimates and the galaxy luminosities of the DOGs have 
been corrected for dust reddening. No reddening correction is given to the other systems. However, those corrections are expected to be 
small. Given the high intrinsic luminosities of the DOG host galaxies, their black holes appear under-massive compared with the z ~ 2 
and local samples. These differences can be explained if the DOGs contain significant young stellar populations which boost the host 
galaxy luminosity. Right: the luminosity expected for the DOGs and other z ~ 2 galaxies after accounting for the fading of their stellar 
populations. An instantaneous burst fading, with Zf ormat ; on at the observed redshift, was applied to the DOGs while passive fading was 
applied to the z ~ 2 galaxies. While the DOGs fade to the local relation, the z ~ 2 systems fade past the loca l relation, suggesting that 
they need to undergo additional stellar mass growth to eventually land on the local relation (Peng ct al. 2006a). Assuming an Eddington 
accretion of ~ 2 M0yr — 1 , the DOGs could easily grow their BHs onto the z ~ 2 relation (after accounting for fading). 



stars in each galaxy ( e.g iKormendv fc Gebhardtl 120011 : 
iFerrarese fc Fo rd 2005), w here t he R-band quantities are 
taken from lBettoni et all (|2003f ) . 

When placing the DOGs into FigureHl we must choose 
what level of dust correction to apply to the BHs and 
host galaxies. For the BHs themselves we show the dust 
corrected lower limits on the masses. For the host galax- 
ies, we use two different prescriptions. In Figure |4] we 
adopt the dust corrections derived from the nuclear re- 
gions. While these are lower limits on the dust extinction 
of the BHs they could be reasonable upper limits on the 
galaxy extinction. Then in Figure [5] we take the other 
extreme, no dust correction on the galaxy host luminosi- 
ties. We discuss both hypotheses below. 

Figure |4] shows that the BH masses measured for the 
DOGs are similar to the z ~ 2 AGNs. However, as a 
function of galaxy luminosity, the DOGs appear in a 
different region of the BH-mass/galaxy-luminosity plot 
compared with the z ~ 2 and local samples. As a func- 
tion of black hole mass, their extinction corrected host 
luminosities are much higher (Figure |4^,). This means 
that either (1) the DOGs have under-massive black holes, 
(2) their host galaxy luminosities are enhanced compared 
with the other galaxy samples, or (3) the dust corrections 
are incorrect. If we take the dust corrections as mea- 



sured, and make the assumption that the DOGs must 
land on the local relation today, they can reach this 
relation by a combination of fading, and/or black hole 
growth. A similar conclusion was reached to explain a 
sample of z ~ 2 SMGs, also with smaller than expected 
black holes for their h ost stellar masses (estimated from 
rest-frame NIR fluxes, IBorvs et al]|2005f) 

The simplest way to explain these differences is if the 
DOGs have larger fractions of very young stellar pop- 
ulations that fade rapidly with time. The DOGs may 
have had significant recent star formation, contributing 
to their high luminosities and large dust reservoirs. In 
contrast, the local galaxies generally have old stellar pop- 
ulations that have already faded from their high-z lumi- 
nosities. If we assume that the DOGs fade as an in- 
stantaneous burst from the time of their observation at 
z ~ 1.5 (e.g. Ivan Dokkum fc Franxl [200l|) . then in fact 
they will fade onto the local BH mass galaxy luminosity 
relation, as shown in Figure^. This fast fading scenario 
is appropriate if mos t of the mass is formed at z ~ 1.5. 

Peng et al.l ()2006aD provide another possible evolution- 
ary scenario for the DOGs on the BH-mass/galaxy- 
luminosity plot. Following the prescription that they 
used for the z ~ 2 AGN in Figures [4] and [5j we can 
assume instantaneous star formation at z = 5, followed 
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Fig. 5. — Left: same as Figure [4] only now no dust correction has been applied to the DOG host galaxy luminosities. Right: now shows 
the DOGs fading passively with an assumed epoch of star formation at redshift, 2 = 5 — the same prescription applied to the other z ~ 2 
AGN hosts. In this scenario, the bulk of the DOGs again fade to the local relation. 



by passive evolution. The essence of this model lies in the 
observation that dusty galaxies with young stellar popu- 
lations generally are more luminous than dust-free galax- 
ies with old populations of equal stellar mass. Therefore, 
passively fading the observed star light at z 2, without 
correcting for extinction, should give an upper limit on 
the light produced by the dominant stellar mass compo- 
nent by z — 0. We show the result of this model in Figure 
where the galaxy stellar luminosities are significantly 
fainter than Figure[4l because we do not correct for dust 
extinction in the galaxy. Even under these assumptions, 
three of the four DOGs still have larger galaxy luminosi- 
ties at a g i yen BH mass compared with the local and 
iPeng et all (|2006al ) high-z samples (Figure [Sh.) . If we 
then apply a passive fading rate, as was applied to the 
other z ~ 2 AGN, the bulk of the DOGs again fade to 
the local relation (Figure [5b). 

In either dust correction scenario, the bulk of the 
DOGs fade onto the local BH- mass/galaxy-luminosity 
relation by today. This behavior is different from the 
z ~ 2 quasar host galaxies which actually fade past the 
local relation, even when the fading applied is only pas- 
sive fading. If the z ~ 2 AGN hosts faded as instanta- 
neous bursts (as has been done with the DOGs in Figure 
3) they would be even further off of the lo cal relation. 
This ph enomenon is discussed extensively in IPeng et al.l 
(2006a); the z ~ 2 AGN host galaxies need to grow in 
stellar mass to land on the BH- mass/galaxy-luminosity 
relation today. Peng et al. estimate that a specific star 
formation rate of 1.2 Gyr -1 would place the quasar hosts 
onto the local relation by z — 1 or, for longer star for- 
mation time scales, the specific star formation could be 
as low as 0.6 Gyr -1 . 



One possible explanation for why the DOGs behave 
differently from the high-z comparison samples is that 
we have under-estimated the BH masses by under- 
estimating the dust corrections on Ha. As a check on the 
BH masses measured from the OSRIS data, we use the 
5 Ks Chandra X-ray observations of the Bootes field to 
make another BH mass estimate. These observations are 
relatively shallow but two of DOGs in our sample are de- 
tected with > 10 counts. For these galaxies, we derived 
the 2 — 10 keV luminosities using an absorbed power- 
law model as the proxy. We fixed the spectral index, 
Gamma= 2, and derived the column density from the 
hardness ratios, taking Galactic columns into account. 
A summary of the fitted parameters is given in Table 
03 We then use the empirical relationship between AGN 
X-ray luminosity and optical continuum flux given by 
(jMaiolino et aljl2007t ). 



log(L5100) 



1 °g( i 2-iokeV) ~ U - 78 



0.721 



(5) 



to derive a second estimate of L5100. Using the X-ray 
derived £5100, we recalculate the BH masses from Equa- 
tion [5] The X-ray derived BH masses are given in Table 
[5j In both cases the BH masses derived from the X- 
ray observations are larger than the lower limits derived 
from the optical OSIRIS data (with large uncertainties). 
The X-ray derived BH mass for DOG2 is estimated to 
be 1.5 — 2.5 times the lower limit from the OSIRIS data, 
while DOG3 has a BH mass about 2 — 3 times larger. 

While the X-ray results push DOGs 2 and 3 to higher 
black hole masses, it is not clear how the estimated host 
galaxy luminosities should change, if at all. The X-ray re- 
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TABLE 5 
X-ray Derived Properties 



Object Name 


X-ray 


0.2-10 keV Flux a 


0.2-10 keV Lum. b 


HI Column 


2-10 keV X-ray Lum. c 


Est. BH Mass 




Counts 


[10 -14 ergs cm -2 s _1 


[10 44 ergs s" 1 ] 


10 21 cm" 2 


[10 44 ergs s" 1 ] 


[10 6 M ] 


DOG2 


16 ±5 


4.4 ± 1.5 


4.2 ± 1.5 


6.1 


2.1 ±0.7 


215 ± 72 


DOG3 


18 ±6 


4.9 ± 1.7 


8.8 ±2.9 


0.6 


3.9 ± 1.3 


3040 ± 850 



Obscrvcd-framc 
^ Rest-frame, uncorrected for internal absorption 
c Rest-frame, corrected for internal absorption 

suits suggest that Ay has been underestimated for these 
two DOGs; not a surprise given the week constraint from 
the TripleSpec spectra on the Ha/H/3 ratios. If we as- 
sume that the additional reddening is primarily affecting 
the central BH, then the AGN (even though they are 
intrinsically brighter than the OSIRIS-based estimate) 
should continue to contribute only a small fraction of 
the rest-frame optical light of the galaxies, because of 
the higher extinction. Thus the host galaxy luminosities 
could remain unchanged. After accounting for luminos- 
ity evolution within the host galaxies (see Figures Hp and 
[5b), the X-ray derived BH masses place the DOGs closer 
to their z ~ 2 counterparts in the BH mass / galaxy lumi- 
nosity plane, although they are also still consistent with 
fading onto the local relation. 

In all of these scenarios the black holes within the 
DOGs may continue to grow. With roughly Eddington 
accretion rates of 2 M0/yr, the DOGs could easily reach 
the z ~ 2 relation in 50 — 100 Myr. Of course if their BHs 
were to grow onto the z ~ 2 relation, the DOGs would 
then also need to continue to grow their stellar mass to 
reach the local relation by today. 

5.2. Star Formation Rates, Metallicities and 
Kinematics 

The apparent (i.e., uncorrected for reddening) star for- 
mation rates of the DOGs in our sample are surprisingly 
small, < 2 M yr -1 . For instance, these rates are a factor 
of 10 — 100 times smaller than the Ha derived star forma- 
tion rates of z ~ 2 SMGs, al so uncorrected for re ddening 
(Figure E SMG data from lSwinbank et al.ll200l . While 
the extinction corrections for DOGs should be large — 
thought to exceed a factor 70 in one of the sample mem- 
bers — even the reddening corrected star formation rates 
are relatively modest. The rates are estimated to be at 
the LIRG level rather than 100 - 1000 M Q yr _1 expected 
in extinction corrected star formation rates of ULIRGs 
and SMGs. 

These low star formation rates may be understandable 
in a scenario where the DOGs are observed in a post 
merger phase and star for mation has begun to shut down. 
iNaravanan et al.l (|2009b| ) attempts to simulate galaxies 
with the unusual selection criteria of D OGs by inducing 
merger s in massive gas rich disks. The INaravanan et al.l 
(2009b) merger simulations actually produce objects that 
would be selected as DOGs. In the simulations, AGN 
dominated DOGs only appear after peak star formation. 
In fact, it is precisely this AGN activity that allows these 
objects to remain DOGs even when the star formation 
rate has dropped below the ULIRG level. 

The observed kinematics of the DOGs can not rule 
out recent merging activity. There is no evidence for 
well ordered rotation in extended star forming disks, and 
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Fig. 6. — The Ha derived SF Rs of the DOGs and a compari- 
son sample of z ~ 2 SMGs from ISwinbank et al.l II2004I) . plotted 
against Ljr. No dust attenuation correction has been applied to 
the SFRs of either sample. Clearly, these four DOGs show signifi- 
cantly less SF than the SMGs despite having similar IR luminosi- 
ties. Note^LfDesWmates of the DOGs are based on 24 fim fluxes 
and Chary & Elbaz (2001) models. They are only meant to be 
representative and do not reflect actual measurements of the FIR 
fluxes of the DOGs. 

the the narrow Ha lines show kinematic offsets from the 
broad Ha lines of 150 — 500 km s _1 . 

It is more difficult to explain the colors and lumi- 
nosities of the DOGs without mergers, but some form 
of rapid gas accret ion from cold flows might work (e.g. 
IBrooks et al.M2009h . However, the metallicity of the star 
forming gas appears to be metal rich with 12 + log(0/H) 
~ 8.6 (solar). Thus, a scenario of pristine gas infall is less 
likely than a scenario of a post merger system where pre- 
vious star formation has polluted the gas with metals. 

6. CONCLUSIONS 

We have obtained high spatial resolution IFU obser- 
vations of four z > 1 extremely dust obscured galaxies 
with the spectroscopy centered on the Ha emission line. 
The Ha emission of all four DOGs are dominated by a 
broad-line component with velocity line widths in excess 
of 2000 km/s. In contrast, the spatially offset narrow Ha 
lines have velocity dispersions < 600 km/s. 

We use the broad Ha line widths and flux measure- 
ments to estimate the BH masses of these systems. After 
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corrections for dust obscuration, the BHs of the DOGs 
are small for their host galaxy luminosity compared with 
z ~ 2 and local unobscured AGN. Much of this offset can 
be explained if the DOGs have significant young stellar 
populations that fade rapidly. Assuming fading from an 
instantaneous burst, the DOGs fade to the local BH- 
mass/host-luminosity relation by today. Alternatively 
the dust correction we have applied to the DOG host 
luminosities could be high. Assuming no dust correc- 
tion for the galaxy luminosity and a passive fading rate, 
the DOGs still evolve onto the local BH-mass/galaxy- 
luminosity relation. 

Interestingly, the z ~ 2 AGN actually fade past the lo- 
cal relation even with only a passive fading applied. Thus 
the comparison samples need to grow in stellar mass to 
reach the local relation today. If the DOGs are in a BH 
growth phase, they could easily reach the z ~ 2 relation 
(after accounting for fading) by growing their BHs at an 
Eddington rate for 50 — 100 Myrs. If they were to do so, 
they too would need to grow in stellar mass to reach the 
local relation. 

The SFRs of the DOGs are also surprising low; af- 
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ter correcting for dust obscuration the SFRs are still 
< 100 M Q yr^ 1 . Thus, without the presence of an AGN 
these galaxies would be classified as LIRGs rather than 
ULIRGs. A merger scenario in which star formation 
peaks prior to the AGN growth phase could explain the 
observations. The measured kinematics of the DOGs can 
not rule out a merger scenario. 
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